Axion isocurvature fluctuations with extremely blue spectrum 
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We construct an axion model for generating isocurvature fluctuations with blue spectrum, niso = 
2 — 4, which is suggested by recent analyses of admixture of adiabatic and isocurvature perturbations 
with independent spectral indices, nad 7^ ^^iBo• The distinctive feature of the model is that the 
spectrum is blue at large scales while scale invariant at small scales. This is naturally realized by 
the dynamics of the Peccei-Quinn scalar field. 



I. INTRODUCTION 

Large-scale structures of the universe, such as galax- 
ies and clusters of galaxies, have formed through grav- 
itational instabilities, initiated by the primordial seed 
density fluctuations, which were created during infla- 
tion. The simplest initial condition seeded these inho- 
mogeneities is the (almost) scale-invariant adiabatic cur- 
vature perturbations. They can fit to very precise mea- 
surements of the cosmic microwave background tempera- 
ture and polarization anisotropics, large-scale structures, 
and supernovae P, Q, Q- K is usually realized by the 
single-field inflation where the inflaton fluctuations are 
responsible for the adiabatic perturbations. 

Generally, there will exist other light fields whose fluc- 
tuations during inflation become isocurvature perturba- 
tions Therefore, the admixture of isocurvature and 
adiabatic fluctuations could be what really happened in 
the early universe. Observational analyses with the as- 
sumption that the spectral indices of the adiabatic and 
isocurvature fluctuations are the same, nad = niso — 1, 
have revealed that the contributions from the isocurva- 
ture perturbation should be small [J, [H| ■ 

However, there is a priori no reason for the isocurva- 
ture fluctuations to have (almost) scale-invariant spec- 
trum. In fact, more general analyses with indepen- 
dent spectral indices of adiabatic and isocurvature modes 
based on recent observations result in the favor of much 
more contribution of the isocurvature component with an 
extremely blue tilt (niso ~ 2 - 4) [i,0,ll- 

In this article, we provide a concrete model for gen- 
erating isocurvature fluctuations with an extremely blue 
spectrum for the first time.^ It is the axion model in 
supersymmetry (SUSY) [l^]. Since the axion is a good 
candidate of the cold dark matter of the universe, and has 
nothing to do with radiation, it gives rise to uncorrelated 
isocurvature fluctuations. The axion isocurvature fluctu- 
ation is usually expressed as 5a/ a ~ H/(2'KFaO), where 
H is the Hubble parameter during inflation. Fa the axion 



decay constant, and 9 a misalignment angle. The key to 
produce the blue spectrum is that we promote Fa as a dy- 
namical field (p = "Fa" , which initially has a large value 
~ Mp, evolves toward smaller values, and stops at ~ Fa 
during inflation. It is realized very simply and naturally 
in the SUSY axion model, and we can obtain extremely 
blue spectrum such as riiso ~ 4 at large scales, which is 
connected to the scale- invariant spectrum at small scales. 

The structure of the article is as follows: In the next 
section, we explain the essence to generate the extremely 
blue spectrum in a simple model, reduced from the con- 
crete model that we provide based on SUSY in Sec. HI. 
We then show the dynamics of the fields which leads to 
the favorable spectrum in Sec. IV. Our conclusions are 
given in Sec.V. 



II. HOW TO GET THE BLUE SPECTRUM 

Let us consider a toy model of a complex scalar field 
<&, whose energy density is negligible during inflation. 
Fluctuations in the phase direction give rise to an isocur- 
vature perturbation, while fluctuations in the radial di- 
rection are negligibly small due to large effective mass in 
that direction as shown shortly. Thus, the isocurvarture 
fluctuation is given by 



H 

2-KLpe' 



(1) 



where we denote $ — ipe^^ I\f2. Since the Hubble pa- 
rameter during inflation is (almost) constant, it is the 
decreasing amplitude of Lp that makes the isocurvature 
perturbation blue tilted. When the field Lp has mass of 
0{H), it can roll down in the potential during inflation, 
and, in addition, its fluctuation Sip is suppressed. The 
reduced potential is given by 



when (p has a large field value, and c ^ 
Then the tp field obeys the equation 



(2) 

0(1) is a constant. 



^ The possibility to obtain isocurvature fluctuations with some de- 
viation from scale-invariant was investigated in Ref. [9|. 



p + 3Hp + cH^p = 0, 



(3) 
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which has a solution of the form ip o< e with 



2 2V 9 ' 



(4) 



for < c < 9/4. Since the isocurvature fluctuation is 
estimated as 
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(5) 



its spectral index is given by 
dhiA^ 



dink =2A = 3-3yi--c. (6) 



Therefore, we obtain the blue spectrum, even extremely 
blue such as niso = 4 for c = 9/4. As shown explicitly 
in the following sections, the field ip eventually settles 
down in the minimum of the potential placed at if Fa- 
Thereafter the isocurvature flucutation becomes scale in- 
variant. 



III. AXION MODEL IN SUSY 

The axion [ll|, [l^l is a Nambu-Goldstone boson asso- 
ciated with the Peccei-Quinn (PQ) symmetry breaking, 
and is the most natural solution to the strong CP prob- 
lem in QCD 13]. The PQ symmetry breaking scale Fa 
is astrophysically and cosmologically constrained within 
the range of lO^" - lO^^ Qgy [IJ]. The axion can be cold 
dark matter for the higher values. 

Let us consider a concrete model of the axion in SUSY. 
We take the following superpotential : 



(7) 



Here $+, and $o are chiral superfields with PQ 

charges +1, —1, and 0, respectively, and h is a constant 
of 0(1). The scalar potential is obtained as 



VsUSY 



F' 



i<f-r)i<foi 



(8) 

where we denote the scalar components with the same 
symbols as the superfields. One can easily see the exis- 
tence of the flat direction which satisfles 



= f: 



$n = 0. 



(9) 



In addition, SUSY breaking effects lift the flat direction 
by soft mass terms 

Kn = rr4|$+|2 + m2_|$_|2+mg|$oP, (10) 



^ One obtains the damping oscillating solution for c > 9/4. Since 
it does not suit for our purpose, we only consider for c < 9/4. 



at low energy scales, where m+, m_, and tuq are of 
O(TeV), as well as the so-called Hubble-induced mass 
terms during inflation. 



(11) 



where c+, c_, and Cq are positive constants of 0(1), which 
stem from the supergravity effects [15i, ilS,]."^ Notice that 
no Hubble-induced A terms will appear due to PQ sym- 
metry. We assume 7? <C in order not to destroy 
the flat direction ([9]). Taking into account the fact that 
4>o = and rrii <^ H (i = ±, 0) during inflation, we only 
consider the potential of the form 



V ^h^ - F^f + c+H^\<S>, 



-c_i/2|$_|2. (12) 



Owing to the Hubble-induced mass terms, the minimum 
of the potential is given by 



1/4 



1/4 



Since it is symmetric between and we consider 
|$+| > |$„| without loss of generality.'* 

Now we must identify the axion field a. Rewriting <^± 

as 



$±.i=^,exp(^e±), 



we can define the fields a and b as 



(^^+^2_)l/2" 

(^^+^2_)l/2«+ (^^+<^2)l/2< 



(14) 

(15) 
(16) 



From Eqs.® or HI]), the potential V{h) for the field h is 
obtained as 



V{h) = -h^Fl^+p^ cos 



(^^ 



2 ^l/2 



(17) 



which implies that the mass of b is given by ~ h{if']^ + 
ip^J)^/"^. Since the field value is ~ Mp initially, and 
decreases until it reaches to Fa, as shown in the next 
section, mf, ^ H during inflation and hence the b field 
quickly settles down into the minimum of the potential. 
On the other hand, the potential for the a field is flat 



^ The coefficients of the Hubble-induced mass terms are deter- 
mined as Ci ~ 3(1 — yi) (i = 0, ±) for the non-renormalizable 
coupling in Kahler potential 5K = yi\^i\^\I\^ /Mp, where / is 
the inflaton and yi's are the coupling constants. 

* Notice that, as shown shortly, the amplitude of the isocurvature 
fluctuation is determined by the larger between |3>+| and |$— |, 
so the spectrum cannot be red tilted. 
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and we can regard it as the axion. During inflation, the 
quantum fluctuations of a develop as 



5a 



(18) 



where iy9_|_ ^ while 56 ~ because it is very massive, 
TOft ^ i/. Thus, (5a_ ~ — (iy9_/iy94-)(5a_(-. Therefore, we 
have 



5e± 



5a± 



± 



H 



(19) 



The crucial point is that the amplitude of the fluctuation 
is determined solely by the larger field value tp^. Also 
notice that the fluctuations in the radial directions Sip-^- 
and Sip^ are both suppressed due to large curvatures in 
their potentials which stem from the flrst term in Vsusy 
[Eq.®] and the Hubble-induced mass terms [Eq. fTT]) ]. 
The axion isocurvature perturbation is given by^ 



Sa 



Tin 



5n~, 



6a 



H 



(20) 



where Ua and denote the number densities of the axion 
and photon, respectively, and we use Eqs. p^ and ^TE\i in 
the last equality. Therefore, the isocurvature fluctuation 
is written as 



Al{k)^A, 



27rV+' 



k—ko 



(21) 

Recent analyses of the admixture of adiabatic and isocur- 
vature perturbations with independent spectral indices, 
n-nd 7^ n-iso, reveal that the isocurvature contribution can 
be as large as the adiabatic mode at the pivot scale fco and 
the blue spectral index of^the isocurvature fluctuation is 
favored such as riiso ~ 4 . 



IV. DYNAMICS OF THE FIELDS AND 
ISOCURVATURE FLUCTUATIONS 

As shown in the previous section, the amplitude of 
the isocurvature fluctuation is solely determined by the 
larger field value iy9_|_ with the constant Hubble param- 
eter during inflation, H ~ const. We therefore need to 
investigate the dynamics of only. It is reasonable to 
consider that the fields slide only along the fiat direc- 
tion, so that if- = 2Fg/Lp+, thus the potential can be 
approximated as 

V c \c+H\l + 2c-H^Ft^ ~ \c+H\l, (22) 



where the last equality holds when tf+ has a large field 
value. ^ Now we must just follow the same argument dis- 
cussed in Sec.IL Since the field obeys the equation 

p+ + iHiP+ + c+H^ip+ = 0, (23) 

whose solution is given by the form cp^ oc e^^^* with 



the isocurvature fluctuation is obtained as 



5a 



H 



so that its spectral index becomes 
dluAE 



cx e 



2XHt 



1 



din fc 



2A = 3-3a/1--c^ 



(24) 



(25) 



(26) 



Therefore, we obtain the blue spectrum with 1 < rtjso < 4 
for < c+ < 9/4. The prominent feature of this model is 
that ip^ eventually settles down to the minimum of the 
potential. 



1/4 



V2(^] Fg, (27) 



and hence we have scale-invariant spectrum afterwards, 
smoothly connected from the blue spectrum at large 
scales. The e-folds during the field evolving from ~ Alp 
to ~ Fn are estimated as 



1, fMp 



(28) 



which gives Nuuc ^ 10 for Fa = lO^^ Q^y ^j^j A = 3/2 
(c+ = 9/4), for example. 

In order to confirm what we have obtained above, we 
solve numerically the equations for and $_ with the 
potential p2| . For the sake of numerical calculations, 
we decompose the field into real and imaginary parts as 
^± — i4>± + 'i4'±)/V2, which leads to the following equa- 
tions: 



2 I aR 



' 2 



21 il 



® For large 95+, kinetic terms are reduced to the normal one as 
^ The actual observable is 5cdm = {^a/f^CDM)Sa oc 9+, where / 4P^\ 1 
the i^xinn densitv na.ra.meter is nc G^. . — > f),, :o^f)^^ ld^ = —[14- — 1 r),,Ln^r)^LD^ ~ — 



the axion density parameter is Ha oc 9^. — y df^ipid^'' ipi = — ( 1 H ^ J i9^i^+9^(/3+ ~ —di_iip-\-d^ ip^ 
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+ + 0«] = 0, (30) 

+ (<^?0'_ + <^i0^) 0i] = 0, (31) 
+ (</'?0'_ + 0^0^) = 0. (32) 
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FIG. 1: Evolution of the fields $+ (upper thick lines) and 
(lower thin lines) for c_ — 9/4 and c+ = 9/4 (niao = 4, red 
solid), 2 (niso = 3, green dashed), and 5/4 (niso = 2, blue 
dotted-dashed). The inset shows the minima where the fields 
settle down. 

Some of the examples are shown in Fig. [TJ Here, we set 
the initial condition as |$+(0)| = Mp. The initial values 
for and the phases d± are taken so as to stay along 
the flat direction ([9]). We also take h — 1 and c_ =9/4. 
In the figure, we show the results for c+ — 9/4, 2, and 
5/4, which correspond to riiso = 4, 3, and 2, respectively. 
One can see that the amplitude of the field (p+ decreases 
exponentially, and eventually stays at the constant value, 
which coincides to Eq. ([771) . This is very attractive, since 
there is no blow-up of the spectrum at smaller scales, 
while having extremely blue tilt even as large as riiso = 4 
at large scales over a few orders of magnitude. Notice 
that the result remain unchanged even if we vary the 
Hubble parameter provided that H ^ Fa; here, we take 



a particular value as H/Fa = 10 ^. 

Finally we comment on the initial amplitude of the 
fields. The ip^ should be at large field values in the be- 
ginning. One of the simple mechanism to realize this sit- 
uation is to consider pre-infiation, where the pre-inflaton 
and the have nonrenomalizable coupling in the Kahler 
potential so as to give rise to a negative Hubble-induced 
mass term during pre-inflation. In this case, the initial 
condition will be ip^ ~ Alp. 

V. CONCLUSIONS 

We have proposed the concrete model for generating 
isocurvature fluctuations with extremely blue spectrum 
for some range of the scale. It is based on the axion 
model in supersymmetry. The supergravity effects raise 
the Hubble-induced mass terms in the potential of the 
ip± fields. These Hubble-induced mass terms play two 
roles. One is that they suppress the fiuctuations in the 
radial directions, Sip±. The other is to make the fields 
evolve during inflation. In particular, the field value of 
(/?+ determines the amplitude of the axion isocurvature 
perturbation: the blue tilt is due to the dynamics of ip+ 
moving from the large initial value (~ Mp) down to the 
PQ symmetry breaking scale Fa during inflation. De- 
pending on the coefficient of the Hubble-induced mass 
term, c+, we can obtain 1 < niso < 4. The prominent 
feature of this model is that the blue spectrum is real- 
ized only while ip^ is evolving and after it settles down 
into the potential minimum the spectrum becomes scale 
invarinat. 

The actual scale L where the spectrum is blue is de- 
termined by e- folds N' after (p+ settles down to Fa- For 
example, niso > 1 at L ^ 1 Mpc for N' ~ 47 assuming 
that the present Hubble radius corresponds to iV ~ 55. 
Observations of large-scale structures, or even PLANCK, 
could see the existence of the isocurvature fluctuations 
with a huge blue tilt, which may approve our model in 
the near future. 
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